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Abstract. We have used the IRAM 30-m and FCRAO 14-m telescopes to observe the molecular clumps associated 
with 12 ultracompact (UC) Hll regions in the J=6-5, 8-7 and 13-12 rotational transitions of methyl-acetylene 
(CH3C2H). Under the assumption of LTE and optically thin emission, we have derived temperature estimates 
ranging from 30 to 56 K. We estimate that the clumps have diameters of 0.2-1.6 pc, H2 densities of 10^-10® cm"'', 
and masses of 10^-210* Mq. We compare these values with those obtained by other authors from different 
molecular tracers and find that the H2 density and the temperature inside the clumps vary respectively like 
riHa oc R"^'^ and T oc i?""'^, with R distance from the centre. We also find that the virial masses of the clumps 
are ~3 times less than those derived from the CII3C2H column densities: we show that a plausible explanation is 
that magnetic fields play an important role to stabilise the clumps, which are on the verge of gravitational collapse. 
Finally, we show that the CH3C2H line width increases for decreasing distance from the clump centre: this effect 
is consistent with infall in the inner regions of the clumps. We conclude that the clumps around UC Hll regions 
are likely to be transient (~10^ yr) entities, remnants of isothermal spheres currently undergoing gravitational 
collapse: the high mass accretion rates (' 
clumps. 



^10 Mq yr ^) lead to massive star formation at the centre of such 
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1. Introduction 

While substantial progress has been done in recent years 
towards a better description of the formation of low-mass 
stars, much remains to be understood of their massive 
counterparts. In fact, several problems hinder the study 
of high-mass star formation: massive stars {M > 10 Mq) 
are more distant, interact more strongly with their envi- 
ronment, and have shorter evolution timescales. After the 
protostellar phase, they reach the ZAMS and create an 
Hii region in the circumstellar medium before the accre- 
tion phase is finished. When the star is very young (i.e. 
< 10^ yr), this ionised region is small (i.e. < 0.1 pc for 
an O star) with respect to the typical diameter of an Hii 
region associated with a more evolved OB star (i.e. ~1- 
10 pc), and is named ultracompact (UC) Hii region. Given 
the strong interaction between early type stars and their 
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surroundings, it is important to have knowledge of the 
physical characteristics of their parental molecular clouds. 
With this in mind, Churchwell et al. (1990), Cesaroni 
et al. (1991), and Hofner et al. (2000) observed respec- 
tively NH3, C^'^S, and C^'^O towards a few UC Hii regions, 
thus assessing the existence of molecular clumps surround- 
ing such UC Hii regions. They estimated masses of 10^- 
10"* Mq, densities of 10'*-10^ cm~'^, and diameters of 0.4- 
1 pc. In another study, Olmi et al. (1993) detected various 
CH3CN rotational lines in a limited sample of UC Hii re- 
gions, which demonstrated the existence of even hotter, 
denser cores inside the more extended clumps observed 
in the C^^'S lines. Notwithstanding these results, a robust 
temperature estimate for the molecular clumps surround- 
ing UC Hii regions was still missing. In fact, on the one 
hand tracers such as NH3 and C^^O provide temperature 
estimates based on 2-3 transitions only and hence prone 
to relatively large uncertainties; on the other, CH3CN is 
an excellent "thermometer" but arises from a region much 
smaller and denser than the 1 pc clumps seen in C^^'S and 
C^^O. Hence, the main goal of the present study was to 



2 



Fontani et al.: CH3C2H in UC Hll regions 



derive a good temperature estimate for the medium den- 
sity molecular surroundings of UC Hii regions. 

The CH3C2H species is especially suitable to our pur- 
poses. This is a symmetric-top molecule (see Townes & 
Schawlow 1975) with a small dipole moment (/x = 0.75 
Debye, DubruUe et al. 1978), which makes thermalisation 
easy even at densities as low as 10'* cm^'^. We thus expect 
CH3C2H to be in "local thermodynamical equilibrium" 
(LTE) conditions in the clumps found by Cesaroni et al. 
(1991) around the UC Hii regions. Moreover, the large 
number of transitions observable in the same bandwidth 
makes the temperature estimate very accurate. For in- 
stance, Bergin et al. (1994) successfully used this method 
to measure the temperature of dense cores in giant molec- 
ular clouds, deriving temperatures of ^ 40 K. 

In Sect. H we describe the observations and data reduc- 
tion procedure, in Sect. |^ we summarise the observational 
results and derive the physical parameters of the clumps, 
in Sect. ^ we discuss the implications of the results ob- 
tained. The conclusions are drawn in Sect. ||. 

2. Observations and data reduction 

2.1. I RAM 30-m telescope 

The data from the IRAM 30-m telescope were obtained 
in the period from August 24 to 29, 1997. The source hst 
is shown in Table 1. The velocity of each source listed 
in Table 1 is taken from previous observations of various 
molecular transitions, and it is only an indicative value of 
the real velocity of the gas from which the CH3C2H lines 
arise. We simultaneously observed the (6-5), (8-7) and 
(13-12) rotational transitions of CH3C2H in the 3 mm, 
2 mm and 1.3 mm bands respectively. The half-power 
beam width (HPBW) of the IRAM 30-m telescope was 
22", 17", and 12" at 3 mm, 2 mm and 1.3 mm, respec- 
tively. 

The alignment between the different receivers was 
checked through continuum cross scans on Jupiter at the 
three frequencies of observation, and found to be accurate 
to within 4". Pointing was checked every 0.5-1 hours by 
means of continuum cross scans on some UC Hll regions 
of the Churchwell et al. (1990) sample. We believe it is 
accurate to within 4". 

The spectrometers used simultaneously are an auto- 
correlator, which covered only the lowest K components 
with high spectral resolution, and a filterbank with low 
spectral resolution, covering all K lines for the (6-5) and 
(8-7) transitions and up to the K=10 component for the 
(13-12) transition. The main parameters of the spectrom- 
eters are listed in Table 2. 

For all sources, the observations were made towards 
the nominal central positions given in Table 1. Only to- 
wards G31.41, we obtained two 3x3 point maps, one with 
12" and the other with 24" spacing. 

We observed using "wobbler switching" , namely a nu- 
tating secondary reflector with a beam-throw of 240" in 
azimuth and a phase duration of 2 sec. The data were cal- 



Table 1. List of the observed sources 



# 


Source 


a(B1950) 


5(B1950) 


wlsr 






(h m s) 


(° ' ") 


(km s '-) 


1 


G5. 89-0.39 


17 57 26.8 


-24 03 56 


10.0 


2 


G9. 52-1-0.19 


18 03 lb. 2 


rjA on AO 

—20 32 03 


4.4 


3 


GlO.47+0.03 


18 05 40.3 


-19 52 21 


67.8 


4 


GlQ.30-0.15 (*) 


18 05 57.7 


-20 06 26 


13.5 


5 


GlO.62-0.38 


18 07 30.7 


-19 56 29 


-3.1 


6 


G19.61-0.23 


18 24 50.3 


-11 58 33 


41.6 


7 


G29.96-0.02 


18 43 27.1 


-02 42 36 


98.0 


8 


G31.41-h0.31 


18 44 59.2 


-01 16 07 


97.0 


9 


G34.26-h0.15 


18 50 46.1 


+01 11 12 


58.0 


10 


G45.47-h0.05 


19 12 04.3 


+11 04 11 


62.0 


11 


W51D (*) 


19 21 22.3 


+14 25 15 


60.0 


12 


IRAS 201264-4104 


20 12 41.0 


+41 04 21 


-3.5 



(*) not observed with the FCRAO telescope 

Table 2. Spectrometers used for the IRAM 30-m obser- 
vations 



Band 


Spectrometer 


Channel Spacing 


Total Band 






(MHz) 


(MHz) 


3 mm 


filterbank 


1.0 


256 




autocorrelator 


0.078 


70 


2 mm 


filterbank 


1.0 


256 




autocorrelator 


0.078 


70 


1.3 mm 


filterbank 


1.0 


512 




autocorrelator 


0.078 


140 



ibrated with the chopper wheel technique (see Kutner & 
Ulich 1981). We checked the calibration by measuring the 
intensities of two planets observed at the same frequen- 
cies as our sources. The main beam brightness temper- 
ature, Tmb, and the flux density, F^, are related by the 
expression F^(Jy) =4.7 Tmb(K). 

2.2. FCRAO Telescope 

Observations of CII3C2H (6-5) were carried out be- 
tween March 16 and May 18 1997 with the Five College 
Radio Astronomy Observatory (FCRAO), located in New 
Salem (Massachusetts, USA). The receiver used was the 
QUARRY focal plane array (Erickson et al. 1992), which 
covered the 3mm band (90-116 GHz) and consisted of 15 
separate beams disposed in a 3 x 5 configuration. The tele- 
scope HPBW at the frequency of the CH3C2H(6-5) line 
is 50". The calibration was carried out using the chopper 
wheel technique. The observations were made in position 
switching, moving the antenna in azimuth by 2'. The main 
beam brightness temperature, Tmb, and the flux density, 
Fi,, are related by the expression F,y{Jy) = 21.5 Tmb(K). 

The spectrometer used was an autocorrelator, with 512 
channels and a spectral resolution of 80 kHz. At FCRAO 
we observed the same sample of UC Hil regions observed 
with the 30-m telescope, except G10.30 and W51D. 
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Fig. 1. Spectra of the CH3C2H (6-5), (8-7), and (12-11) rotational transitions towards sources G10.62, G10.47 and 
G45.47. On top, the FCRAO (6-5) spectra are shown 



2.3. Data reduction and fitting procedure 

CH3C2H is a symmetric-top molecule, with rotational lev- 
els described by two quantum numbers: J is associated 
with the total angular momentum, whereas K is associ- 
ated with its projection on the symmetry axis. For each 
value of J, the selection rules A J = ±1 and AK — make 
possible only rotational radiative transitions J + 1 ^ J 
with K < J] consequently, for each J + 1 — > J transition, 
only J + 1 X-lines can be seen. In our observations we 
could detect only lines up to a maximum of K — 4. 

The data reduction made use of the CLASS program of 
the GAG software developed at IRAM and Observatoire 
de Grenoble. In order to derive the line parameters of 
CH3C2H we assumed that all the K components of each 
J -|- 1 — > J transition arise from the same region, and 
are hence characterised by the same peak velocity and 
line width. Then, Gaussian fits to each J -I- 1 — > J spec- 
trum were made with line separations between different K 
components fixed to the laboratory values and line widths 
assumed to be identical. 

3. Results 

With the IRAM 30-m telescope we detected CH3C2H (6- 
5), (8-7) and (13-12) towards all of the 12 sources of our 
sample, whereas with the FCRAO we detected CH3C2H 
(6-5) towards all of the sources but G19.61. 

In Tables 3, 4, 5 and 6 we list the line parameters 
obtained with the data reduction procedure described in 



Sect. 2.3. In all cases, the difference between the line pro- 
file and the Gaussian fit was within the noise. For each 
source we list: the LSR velocity, Vlsr,; the line full width at 
half maximum (FWHM), Avi/2; the emission integrated 
under each line, J Tmb^v; and the la rms noise of the 
spectrum. For Gaussian profiles, the peak temperatures 
are related to J Tmb^v by the relationship: 



7^ max 
-'MB 



Avi/2^ -^^^ 



(1) 



The errors quoted do not take into account the calibration 
uncertainties, but represent only the formal errors of the 
fit procedure. The calibration errors will be discussed in 



Sect. 3.2 



3.1. Size 

In order to assess that the CH3C2H lines originate from 
the 1 pc clumps known to exist around the UC Hii regions, 
we must estimate the size of the CII3C2II emitting re- 
gion. The maps obtained with the QUARRY multi-beam 
receiver on the FCRAO telescope are unsuitable to this 
purpose, because emission was detected only towards the 
central position. Consequently, only an upper limit to the 
source size could be obtained. However, another estimate 
of the size can be derived from the comparison between the 
CH3C2H(6-5) line intensities measured with the IRAM 
30-m and FCRAO telescopes. Assuming Gaussian inten- 
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Table 3. Parameters of the CH3C2H(6-5) lines observed with the IRAM 30-m telescope. The results given here are 
based on the autocorrelator (high resolution) spectra. The errors are given by the fit procedure, and do not consider 
the calibration error on Tmb- ^^1/2 represents the full width at half maximum (FWHM) of the line, a is the noise 
level in K 



Source 




Alli/2 




/■ TMBd^; (K km s 






(7 




(km s"^) 


(km s^^) 


A' = 


1 


2 


3 


4 


(K) 


G5.89 


8.54±0.08 


3.45±0.1 


15.15±0.1 


14.07±0.08 


7.24±0.1 


5.35±0.1 


0.85±0.1 


0.1 


G9.62 


4.41±0.02 


3.87±0.04 


7.45±0.07 


6.23±0.06 


2.82±0.06 


1.95±0.06 


0.23±0.06 


0.04 


G10.47 


66.66±0.02 


6.46±0.07 


6.77±0.1 


5.8±0.1 


2.94±0.1 


2.4±0.1 


0.24±0.1 


0.1 


G10.30 


13.06±0.05 


3.87±0.03 


2.18±0.05 


1.78±0.05 


0.81±0.05 


0.38±0.05 


0.08±0.05 


0.04 


G10.62 


-2.82±0.03 


6.13±0.06 


9.69±0.09 


8.41±0.08 


4.54±0.08 


3.28±0.08 


0.38±0.07 


0.06 


G19.61 


41.74±0.09 


6.32±0.06 


3.42±0.09 


2.93±0.09 


1.48±0.08 


1.10±0.07 


0.1±0.07 


0.05 


G29.96 


97.70±0.01 


2.75±0.03 


1.57±0.05 


1.40±0.04 


0.74±0.03 


0.59±0.04 


0.05±0.04 


0.05 


G31.41 


98.25±0.05 


4.31±0.05 


5.1±0.1 


3.8±0.1 


2.1±0.1 


1.6±0.1 


0.15±0.1 


0.1 


G34.26 


57.59±0.04 


5.38±0.06 


7.05±0.09 


6.64±0.08 


3.44±0.08 


2.18±0.08 


0.10±0.08 


0.09 


G45.47 


62.37±0.09 


3.87±0.05 


3.0±0.1 


2.0±0.1 


0.9±0.1 


0.75±0.1 


< 0.6 


0.1 


W51D 


61.06±0.04 


4.30±0.06 


6.7±0.1 


5.7±0.1 


3.25±0.1 


2.8±0.1 


0.3±0.1 


0.1 


IRAS 20126 


-3.70±0.07 


2.18±0.02 


1.84±0.04 


1.52±0.04 


0.65±0.05 


0.39±0.04 


< 0.2 


0.04 



Table 4. Parameters of the CH3C2H(6-5) lines observed with the FCRAO telescope. G10.30 and W51D were not 
observed 



Source 


ULSR 






I'Tyoidv (K km s"') 




O" 




(km s^^) 


(km s-i) 


K = Q 


1 


2 


3 


(K) 


G5.89 


8.79±0.05 


3.46±0.03 


6.7±0.2 


5.6±0.2 


3.0±0.2 


2.0±0.2 


0.2 


G9.62 


4.2±0.1 


4.2±0.1 


2.9±0.1 


2.6±0.1 


1.2±0.1 


0.6±0.1 


0.1 


G10.47 


66.9±0.2 


5.4±0.2 


2.2±0.2 


1.9±0.1 


0.9±0.1 


0.7±0.1 


0.09 


G10.62 


-2.9±0.1 


5.3±0.1 


4.3±0.2 


3.2±0.2 


1.5±0.2 


1.1±0.2 


0.1 


G19.61(*) 






< 3.0 








0.15 


G29.96 


97.8±0.2 


3.6±0.2 


1.05±0.09 


0.83±0.08 


0.29±0.08 


0.43±0.08 


0.1 


G31.41 


97.66±0.09 


4.2±0.1 


2.2±0.1 


1.9±0.1 


l.OiO.l 


0.8±0.1 


0.09 


G34.26 


57.52±0.04 


4.51±0.07 


3.79±0.07 


3.5±0.1 


1.5±0.1 


l.OiO.l 


0.1 


G45.47 


62.3±0.2 


4.7±0.3 


1.39±0.09 


1.12±0.09 


0.41±0.08 


< 0.3 


0.07 


IRAS 20126 


-3.68±0.06 


1.7±0.07 


0.58±0.03 


0.54±0.03 


0.22±0.02 


< 0.2 


0.05 



(*) undetected source: the upper limits on JTmb have been obtained from a assuming a line FWHM equal to that measured 
with the IRAM 30-m telescope. 



sity distributions for the CH3C2H emitting clumps, one 
finds that the main beam brightness temperatures, Tmb, 
measured with the two telescopes are related by the ex- 
pression 



pfcrao 

-MB 



(9 



2 

fcrao 



A, 



-MB 



e 



(2) 



where is the FWHM of the source, and 830m and Ofcrao 
are the instrumental HPBWs. 

From this equation one obtains the expression for 0s : 



02 

q2 _ fcrao 



C\2 /7^30m /7^fcrao'\ 



/T^SOm /T^fciao\ i 



(3) 



The values of Qs are listed in Table 7 and the corre- 
sponding linear diameters, Z?, are given in Table |. The er- 
rors have been computed from Eq. (|^) with the usual prop- 
agation of the statistical errors and hence depend upon the 
which will be discussed in Sect. 3^. The 
30" and D ~ 0.9 pc. We stress 
that in all cases Qs is consistent with the upper limits ob- 
tained from the FCRAO maps. Also, for G31.41 the value 



errors on Tmb, 
mean values are Q 



derived in this way can be compared with the direct es- 
timate obtained from the emission map made with the 
IRAM 30-m telescope. Figures || and || show the maps of 
G31.41 obtained by integrating respectively the (6-5) and 
(8-7) emission under the K—Q and 1 lines. For a Gaussian 
source, one can compute the deconvolved clump diame- 
ter from the observed FWHM. We find 63 = 29" from 
the GH3C2H(6-5) map and 9^ = 26" from CH3C2H(8- 
7) . These estimates are consistent within the uncertainties 
with that obtained from Eq. (|), 63 = 32" ±8", thus con- 
firming the reliability of the method. No size estimate can 
be derived from the GH3C2H( 13-12) map because emis- 
sion was detected only towards the central position: how- 
ever, the upper limits obtained at offsets of ±12" from the 
centre are consistent with the (13-12) line intensity mea- 
sured towards the (0,0) position and the source diameter 
estimated above. 

We can also compare our results with those derived 
from G^^O (Hofner et al. 2000) and G^^S (Gesaroni et al. 
1991): in the former, the angular diameters are slightly 
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Table 5. Parameters of the CH3C2H(8-7) lines observed with the IRAM 30-m telescope 



Source 








/■ Tyisdv (K km s 






a 




(km s-i) 


(km s"^) 


K = 


1 


2 


3 


4 


(K) 


G5.89 


8.52±0.08 


3.50±0.05 


26.6±0.2 


25.3±0.2 


14.8±0.2 


12.2±0.2 


2.4±0.2 


0.2 


G9.62 


4.67±0.02 


3.67±0.06 


10.08±0.09 


8.34±0.09 


4.20±0.09 


2.99±0.08 


0.80±0.08 


0.07 


G10.47 


66.12±0.03 


7.5±0.1 


18.8±0.2 


11.7±0.2 


6.5±0.2 


6.2±0.2 


3.2±0.2 


0.2 


G10.30 


13.03±0.05 


3.85±0.03 


3.45±0.07 


2.80±0.07 


1.45±0.07 


0.90±0.07 


0.1±0.06 


0.08 


G10.62 


-2.92±0.03 


5.88±0.08 


14.2±0.1 


12.5±0.1 


7.7±0.1 


6.4±0.1 


1.4±0.1 


0.1 


G19.61 


41.19±0.08 


6.62±0.09 


6.2±0.1 


4.7±0.1 


2.8±0.1 


2.5±0.1 


l.OiO.l 


0.1 


G29.96 


97.69±0.01 


3.2±0.1 


3.45±0.07 


2.82±0.07 


1.52±0.07 


1.28±0.06 


0.18±0.06 


0.08 


G31.41 


97.74±0.03 


5.60±0.08 


10.50±0.2 


6.0±0.1 


3.8±0.1 


3.9±0.1 


l.lzbO.l 


0.2 


G34.26 


57.67±0.04 


6.07±0.09 


16.0±0.1 


ll.lzbO.l 


6.1±0.1 


5.3±0.1 


l.lzbO.l 


0.2 


G45.47 


62.49±0.09 


3.5±0.1 


4.0±0.2 


3.5±0.2 


1.8±0.2 


0.9±0.2 


< 0.7 


0.2 


W51D 


60.84±0.04 


5.22±0.1 


15.0±0.2 


10.7±0.2 


7.3±0.2 


6.3±0.2 


1.7±0.2 


0.2 


IRAS 20126 


-3.77±0.07 


2.37±0.06 


3.025±0.07 


2.60±0.06 


1.10±0.06 


0.60±0.06 


0.06±0.05 


0.07 



Table 6. Parameters of the CH3C2H(13-12) lines observed with IRAM 30-m telescope 



Source 


VLSR 


A«l/2 




/■ TMsdv (K km s" 


') 




a 




(km s"^) 


(km s"^) 




1 


2 


3 


4 


(K) 


G5.89 


8.56±0.06 


3.92±0.06 


29.0±0.3 


28.6±0.3 


18.6±0.3 


18.0±0.3 


4.0±0.3 


0.4 


G9.62 


4.78±0.05 


4.16±0.06 


6.7±0.2 


6.0±0.2 


4.2±0.2 


3.6±0.2 


0.5±0.2 


0.2 


G10.47 


66.55±0.1 


8.0±0.1 


15.0±0.4 


16.0±0.2 


9.6±0.2 


9.5±0.1 


3.1±0.1 


0.4 


G10.30 


12.79±0.08 


3.47±0.09 


2.6±0.1 


2.4±0.1 


1.4±0.1 


l.OiO.l 


0.2±0.1 


0.2 


G10.62 


-2.74±0.05 


6.1±0.1 


14.6±0.3 


14.5±0.3 


9.6±0.3 


9.6±0.2 


2.2±0.2 


0.3 


G19.61 


40.8±0.1 


9.1±0.2 


5.7±0.3 


10.0±0.3 


6.5±0.2 


5.2±0.2 


0.4±0.2 


0.2 


G29.96 


98.01±0.05 


3.53±0.07 


3.1±0.1 


2.7±0.1 


1.7±0.1 


1.7±0.1 


0.2±0.1 


0.2 


G31.41 


97.96±0.1 


5.7±0.1 


10.62±0.1 


9.7±0.1 


9.5±0.1 


7.3±0.1 


2.1±0.1 


0.3 


G34.26 


58.28±0.05 


6.79±0.07 


12.1±0.3 


10.88±0.3 


7.9±0.2 


6.9±0.2 


1.8±0.2 


0.2 


G45.47 


62.9±0.2 


3.6±0.2 


2.5±0.3 


3.2±0.3 


2.2±0.3 


1.7±0.3 


< 1-1 


0.4 


W51D 


60.59±0.06 


4.99±0.07 


16.0±0.2 


15.9±0.2 


lO.OiO.l 


9.8±0.3 


1.6±0.3 


0.4 


IRAS 20126 


-3.44±0.1 


3.26±0.1 


1.9±0.1 


1.8±0.1 


0.8±0.1 


0.6±0.1 


O.liO.l 


0.2 



greater (^20%) than those estimated by us, with the only 
exception of G10.47, which is two times larger in C^^O; 
in the latter tracer, instead, the angular sizes are ~50% 
smaller with respect to CH3C2H. We shall come back to 
this point in Sect. 4.2. 



3.2. Temperature and column density 

From the angular diameters obtained from Eq. (^ , we de- 
rived the rotation temperature and total column density 
of the CH3C2H molecules by means of the population dia- 
gram method (see e.g. Hollis 1982, Olmi et al. 1993), which 
assumes the gas to be in local thermodynamical equilib- 
rium (LTE). Such an assumption is believed to work very 
well for CH3C2H (see Bergin et al. 1994) due to its low 
dipole moment. Under the additional assumption of opti- 
cally thin emission, one can compute the beam averaged 
column density Ni of the upper level i from the line in- 
tensity: 



(4) 



Ni _ J TMBdw 

where k is the Boltzmann constant, rji, the beam filling 
factor, V the line rest frequency, 5* the line strength, 



the statistical weight, and fi = 0.75 Debye (DubruUe et 
al. 1978) the dipole moment of the CH3C2H molecule. 

In Table we give the rotational temperatures {T^ot) 
and total CH3C2H source column densities (A'^tot): these 
have been obtained from a least square fit to the data, as 
shown in the "Boltzmann plots" in Fig. |^. We note that 
each column density A^^ in Fig. ^ includes the correction 
for angular resolution effects, assuming gaussian profiles 
for both the instrumental beam and the source. Therefore, 
A^tot corresponds to a mean value over the clumps. The 
temperatures are distributed between 30 and 50 K and 
the column densities are of order 10^'*-10^^ cm~^. The 
goodness of the fits in Fig. ^ confirms that the gas is very 
close to LTE, which in turn implies that the rotational 
temperature derived in this way is an excellent estimate 
of the kinetic temperature, Tk, of the molecular gas. 

The errors quoted in Table |^ take into account calibra- 
tion uncertainties on Tmb of 10% for the 3 mm lines and 
20% for the others. For Trot and A^tot, the uncertainties 
have been estimated allowing for a variation of the values 
of A^i by 10% simultaneously for all of the CH3C2H(6-5) 
K fines and by 20% for the (8-7) and (13-12) fines: this 
takes into account the fact that each sample of J -I- 1 —> J 
transitions is affected by the same calibration uncertainty. 
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-20 - 



20 -20 

ha (arcsec) 

Fig. 2. Map obtained using the IRAM 30-m telescope to- 
wards G31.41 by integrating the CH3C2H(6-5) emission 
under the K=Q and K=\ hnes. Contour levels range from 
1 K km s^^ (corresponding to 3cr) to 7 K km s^^ in steps 
of 2 K km s"'^. The thick line corresponds to the half max- 
imum power contour 



20 




-20 



20 -20 

Acx forcsBcj 

Fig. 3. Same as Fig. |for the CH3C2H(8-7) K=0 and 1 
lines. Contour levels range from 2 K km s~^ (correspond- 
ing to Scr) to 18 K km s~^ in steps of 4 K km s~^ 



Then the errors have been taken equal to the largest dif- 
ference between the values of Trot and A'tot thus obtained 
and the nominal ones. 

As previously explained, the population diagram 
method uses the assumptions of LTE and optically thin 
emission: let us now discuss these in some better detail. 

The LTE hypothesis is acceptable when the H2 density 
in the clump is greater than the critical density (see e.g. 
Spitzer 1978). The maximum critical density for the tran- 




50 100 150 200 



50 100 150 200 



(K) 

Fig. 4. Boltzmann plots for all of the observed sources. 
Filled circles represent the CH3C2H (6-5) data, triangles 
the (8-7), and empty circles the (13-12). The column den- 
sities have been obtained from the IRAM 30-m spectra. 
The straight lines represent least square fits to the data. 
The error bars correspond to the formal errors of the fit- 
ting procedure and do not consider the calibration error 
on Tmb 



sitions observed by us is ^ 10^ cm~^: we shall see in the 
next section that the molecular hydrogen densities turn 
out to be greater than this value, which supports our as- 
sumption of LTE (and hence T.^t ^ Tk) for the CH3C2H 
molecules. 

As for the optically thin assumption, it is worth noting 
that the observed bandwidth covers most of the K com- 
ponents of the CH3I3CCH (6-5), (8-7), and (13-12) tran- 
sitions. One can thus estimate the optical depth from the 
intensity ratio between the lines of the main species and 
those of the isotopomer. However, none of the CH3^^CCH 
lines is detected by us, so that only an upper limit on 
the optical depth can be obtained: assuming an isotopic 
ratio i2c/i3c=50 (see Wilson & Rood 1994), the mini- 
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Table 7. Angular diameters, rotation temperatures, and 
source averaged column densities of the CH3C2H clumps. 
The errors are calculated assuming a calibration error of 
~ 10% at 3 mm and ^ 20% at 2 mm and 1.3 mm (see 
text) 



Source 




T^ot 


Atot 




(arcsec) 


(K) 


(10"cm-2) 


G5.89 


31±8 


48±10 


20±1 


G9.62 


28±7 


34±6 


13±1 


G10.47 


22±5 


46±8 


27±1 


G10.30(*) 


25 


34±7 


5±1 


G10.62 


27±8 


46±10 


21±1 


G19.61(*) 


30 


56±7 


11±3 


G29.96 


45±9 


48±8 


4±1 


G31.41 


32±8 


53±8 


15±2 




29±5(**) 


51±7 


15±3 




26±5(***) 


48±7 


16±3 


G34.26 


40±9 


46±9 


16±2 


G45.47 


35±8 


39±7 


5.0±0.3 


W51D(*) 


23 


47±10 


15±3 


IRAS 20126 


22±7 


30±5 


3.9±0.2 



(*) Sources for which the angular diameter could not be com- 
puted. For G10.30 and W51D Trot and Atot have been com- 
puted assuming a linear diameter equal to the mean value 
(0.9 pc) obtained from the other sources. For G19.61, instead, 
we have assumed the diameter of 30" derived by Hofner et al. 
(2000) from the C"0 emission 

(**) Aiigular diameter measured from the CH3C2H(6-5) map 
of Fig. 

(***) Augular diameter measured from the CH3C2H(8-7) map 

of Fig. r 



mum value turns out to be t=8, which is insufficient to 
demonstrate that the CH3C2H lines are optically thin. 
Nevertheless, the excellent agreement between the linear 
fit and the data in the Boltzmann plots reinforces the as- 
sumption of optically thin line emission: in fact, large opti- 
cal depths are expected to affect mostly the low excitation 
transitions and hence cause a flattening of the plots at low 
energies, which is not seen in Fig. |4[ 

We will examine the optical depth issue in some better 



detail in Sect. 3.4.1 



Another topic which may be discussed in relation to 
the fits in Fig. ^ is the effect of temperature and den- 
sity g radien ts o n the rotation diagrams: we delay this to 
Sects. I3.4.2I and b.4.3. 



3.3. Mass and density 

By means of the parameters derived above, we can esti- 
mate the mass of the clumps and the corresponding virial 
mass needed for equilibrium. For both estimates, spherical 
symmetry is assumed. 



The mass of the clump, Mcd, can be obtained from 
the column density in Table |^: for a spherical, isothermal 
clump one has 



CD 



4 



(5) 



ArcH3C2H 

where d is the distance to the source, ttihs the mass of 
the H2 molecule, ArcH3C2H the abundance of CH3C2H rel- 
ative to H2 , and A'tot is the mean column density over the 
clump, as previously explained. 

The most uncertain value in Eq. (H) is ^ch3C2H- 
Rather than assuming a priori an arbitrary value for a 
parameter which is known to vary by a large amount in 
the interstellar medium, we have decided to obtain an es- 
timate of it by equating Eq. (^) to the mass, Afcont, de- 
rived from continuum measurements of our sources (the 
latter is given in Table ^). To this purpose, we used the 
(sub)millimeter continuum emission observations made by 
various authors (Chini et al. 1986; Cesaroni et al. 1999; 
Hatchell et al. 2000; Hunter et al. 2000) towards most of 
our sources. Following Mezger et al. (1990), the mass was 
derived using these continuum measurements and the dust 
temperatures quoted by the corresponding authors; when 
the latter were not available, we used Trot from Table 0. 
The values of ArcH3C2H obtained by us range from 0.9 10~^ 
to 5.7 10-^ (see Table |) with a mean value of ~ 2.2 10"^. 
We have hence assumed XCH3C2H = 2 10^^, with an mi- 
certainty of a factor ~2.5 and computed the A/cd from 
Eq. (|^) using this value for all sources. The results are 
given in Table ^. 

It is worth pointing out that the continuum measure- 
ments do not necessarily refer to the same region seen in 
the CH3C2H lines: however, in a few cases (G5.89, G9.62, 
G10.47, G31.41) maps in the (sub)millimeter continuum 
are available, so that we could verify the coincidence be- 
tween the CH3C2H and dust emitting regions (~30"). For 
the remaining sources, the continuum observations were 
performed with an instrumental HPBW greater than ours, 
so that the mass estimates in Table ^ are to be taken as up- 
per limits to the mass of the clumps traced by CH3C2H, 
which in turn imply lower limits on the CII3C2H abun- 
dances; however, the relatively small range of variation 
of >'^CH3C2H in Table ^ suggests that our assumption of 
XCH3C2H = 2 10"^ is likely to be a reliable approxima- 
tion. 

The virial mass for a homogeneous spherical clump, 
neglecting contributions from magnetic field and surface 
pressure, can be computed from the expression (see e.g. 
MacLaren et al. 1988) 

Mvir(MQ) = 0.509 d(kpc) es(arcsec) AD^/2(km/s) (6) 

where Aui/2 is the observed line FWHM. In our case we 
took Avi/2 from the CIl3C2H(6-5) lines, for the sake of 
consistency with the derivation of 6s which has also been 
obtained from the CH3C2H(6-5) emission. The values of 
Mvir are given in Table ^ 

In the same table we also give the H2 volume densi- 
ties estimated from Mcd and the clump diameters, D. 
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The uncertainties on these values have been calculated by 
means of the propagation of statistical errors, but they do 
not take into account distance uncertainties. All masses 
are equal to a few ^ IQ'^AIq, with the sole exception of 
IRAS 20126, for which the mass is < 10^ M©. The H2 vol- 
ume densities range from 10^ to 5 10^ cm~'^, greater than 
the maximum critical density of the transitions observed 



(see the discussion in Sect. 3.1) 



3.4. Effects of inhomogeneity and optical depth 

In view of the findings that will be illustrated in Sect. ^, 
before proceeding further it is important to discuss the 
relevance of optical depth effects, temperature gradients, 
and density gradients on the mass (and hence density) 
estimates obtained so far. 



3.4.1. Optical depth 



As explained in Sect. 3.2, we believe that the bulk of the 
observed CH3C2II emission is optically thin. However, one 
cannot exclude the existence of an optically thick region 
close to the centre of the clump, where the density is 
presumably higher than average. Indeed, dense molecu- 
lar cores with diameters of ~0.1 pc are known to exist in 
most of our clumps. Hence the question is how much of 
the clump mass is contained in the optically thick region. 
One can calculate the size of the region over which r > 1 
for the CH3C2H hues. For A^tot = 2 lO^^ cm-^ T = 40 K, 
Avi^2 = 5 km s~^, and a density profile n oc one finds 
that the optical depth of e.g. the CH3C2H(6-5) K=0 tran- 
sition along the line of sight passing through the centre of 
the clump is equal to r ~ Q.1{Rq/R — 1), where Rq is the 
radius of the clump and R the distance from the centre at 
which the optical depth r is achieved. By posing t = 1 one 
finds the radius of the optically thick core, R/Rq — 10. In 
conclusion, the CH3C2H gas becomes optically thick at a 
distance from the centre equal to ~10% of the clump ra- 
dius: this contains ~10% of the total mass. We conclude 
that optical depth effects are negligible for the estimate of 

McD. 

Incidentally, it is worth noting that in case of large op- 
tical depths the values of A'tot and Mqd are to be taken 
as lower limits. On the contrary, Mvir becomes an upper 
limit, because for r ^ 1 line broadening occurs, which 
increases the value of Avi/2 used in Eq. (^. As a con- 
sequence, the ratio Mqd /Mvir which one can derive from 
Table || has to be regarded as a lower limit. This fact re- 
inforces the conclusion attained later in Sect. 4.1. 



3.4.2. Temperature gradient 

The hypotheses under which the method of the popu- 
lation diagrams in Sect. 3^ is applied are that the gas 
is optically thin, in LTE, and isothermal. The first two 
assumptions have already been discussed; as for the lat- 
ter, Mauersberger et al. (1988) have treated the case of 



CP 




RAS 20126 
III 




3 4 5 3 4 5 

ln[Ei(K)] 

Fig. 5. Same as Fig. ^ but with logarithmic scale on the 
abscissa 

a spherically symmetric cloud with temperature and den- 
sity having power-law dependence on the distance from 
the centre. These authors demonstrate that in this case a 
linear correlation should hold between In Ni and In Ei , un- 
like the isothermal case where the correlation is expected 
with Ei. We stress that the different dependence on Ei is 
caused only by the temperature gradient: in fact, for an 
isothermal cloud the usual linear relation between In Ni 
and Ei is recovered, no matter what the density gradient 
is. 

Following the example of Mauersberger et al. (1988), 
in Fig. ^ we present "modified" population diagrams for 
our sources, plotting In Ni versus \nEi. Clearly, the linear 
fits to the points are much less satisfactory than those in 
Fig. ^. We conclude that the CH3C2H emission arises from 
a region in which the temperature profile must be very 
shallow, hence making r=const. a suitable approximation. 

3.4.3. Density gradient 

For an isothermal and optically thin cloud, the source av- 
eraged column density cannot depend on the way the gas 
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Table 8. Mass and CH3C2H abundance estimates from 
the continuum millimeter emission 



10 - 



Source 


m — 

(kpc) 


Mcont 


-^CH3C2H 

(xlO-^) 


Ref. 


G5.89 


4.0 


1500 


5.5 


1 


G9.62 


5.7 


1600 


5.7 


1 


G10.47 


5.8 


12700 


0.9 


2 


G10.30 


6.0 








G10.62 


6.0 


8700 


1.7 


3 


G19.61 


3.5 


2200 


1.4 


3 


G29.96 


7.4 


11200 


1.0 




G31.41 


7.9 


7300 


3.5 


2 


G34.26 


4.0 


9700 


1.1 


3 


G45.47 


8.3 


12900 


0.9 


3 


W51D 


8.0 








IRAS 20126 


1.7 


110 


1.3 


4 



Reference codes: l=Hunter et al. (2000); 2=HatchelI et al. 
(2000); 3=Cliini et al. (1987); 4=Cesaroni et al. (1999) 

distances are from Wood & Churchwell (1989a) for G10.30 
and from Kurtz et al. (2000) for the other sources 



is distributed inside the cloud: hence, density gradients 
leave the estimates of A'tot and Mqd unaffected. On the 
other hand, the virial mass depends on the density pro- 
file, as demonstrated e.g. by MacLaren et al. (1988). For a 
power-law density distribution of the type riHj oc R^, the 
virial mass obtained from Eq. (|^) must be multiplied by 
the factor which is <1 for p < 0. Therefore, the 

values of M^ir of Table ^ derived for a homogeneous cloud 
are to be taken as upper limits. 

In conclusion, the ratio Mcd /Afvir derived by us repre- 
sents an underestimate if the constant density assumption 
is released in the computation of M^ir- This result is im- 



portant in view of the discussion in Sect. 4.1 



4. Discussion 

In the previous section, we have demonstrated that the 
values of the physical parameters derived by us are robust. 
In particular, we have shown that: (i) optical depth effects 
can be neglected for the CII3C2H lines; (ii) temperature 
gradients, if present, must be very "flat" in the CII3C2II 
emitting region; and (iii) density gradients have no effect 
on the estimate of Mcd and the corresponding H2 density. 
In the following, we shall use these physical parameters to 
analyse the stability and structure of the clumps. 



4.1. Stability of the clumps 

From Table ^ one can see that the mean ratio between 
Mcd and Mvir is '-~^3.0±1.7, being ~I only in the cases 
of G19.61 and IRAS 20126. Ahhough such a ratio is only 
marginally >1, it is such for all of the sources but two of 
them: this result seems too systematic to be due only to 
random errors on the quantities of interest. The question 




Fig. 6. Plot of the ratio between the mass estimated from 
the gas column density and that derived assuming virial 
equilibrium for the clumps observed in this study. Full cir- 
cles indicate values computed from CH3C2H (this study), 
empty triangles those derived from C^^'S (Cesaroni et 
al. 1991), and filled squares those obtained from C^^O 
(Hofner et al. 2000). The straight line corresponds to 
Mcd = Mvir 



is whether the ratio AfcD/Mvir can be reduced to unity 
by any means. Let us examine a few possibilities. 

opti- 



As already discussed in Sects. 3.4.1 and 3.4.3, 
cal depth effects and density gradients can only increase 
the ratio Afco/Mvir, as they make A/cd bigger and Mvir 
smaller. 

Another possibility is that the CII3C2H abundance 
is underestimated, which would cause an overestimate 
of Mcd- In principle this cannot be excluded, given the 
large variations of molecular abundances in different ob- 
jects (see e.g. Irvine et al. 1987); however, our estimate 
of XCH3C2H (see Sect. [3^ ) is "direct", as it makes use of 
the continuum emission from the regions of interest and 
it is hence less prone to errors which instead affect values 
based on chemical models. Further evidence for Mcd be- 
ing greater than Mvir comes from the studies of Cesaroni 
et al. (1991) and Hofner et al. (2000) who mapped some 
of our clumps respectively in the C^^'S and C^^O hues: as 
one can see in Fig. |[ in all sources the ratio Mcd /Mvir 
is basically >1, independently of the tracer. It is worth 
stressing that the values of Mcd obtained by Hofner et 
al. (2000) from C^^O are less affected by uncertainties on 
the abundance, which is much better established for C^^O 
than for CH3C2H. 

In conclusion, we favour the hypothesis that the dis- 
crepancy between Mcd and Mvir is real. In this case 
Mcd /Mvir > 1 seems to indicate that the clumps are un- 
stable against gravitational collapse. 

If thermal pressure and turbulence are the only means 
of support against gravitational collapse, indeed Mcd > 
Mvir implies collapse of the clumps on a free-fall time: for 
densities of order 10^ cm~'^ the latter is equal to ~10'' yr, 
an order of magnitude less than the estimated lifetime 
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Fig. 7. Plot of the mass accretion rate versus the source 
bolometric luminosity estimated from the IRAS fluxes. 
Circles and squares indicate estimates obtained respec- 
tively from CH3C2H (this study) and C"0 (Hofner et al. 
2000). Empty circles indicate the sources for which it was 
not possible to estimate the size of the CH3C2H clump (see 
Table 0) . The straight line is the best fit assuming propor- 
tionality between mass accretion rate and luminosity and 
corresponds to Macc(-W0yr~^) = 2.1 lO~^LFiR(i0) 



of the UC Hii regions embedded in the clumps (Wood 
& Churchwell 1989b). This suggests that magnetic fields 
might play a role in stabilising the cloud against gravita- 
tional collapse. From the virial theorem one can estimate 
the strength of the magnetic field, B, required to this pur- 
pose: following McKee et al. (1993) and neglecting the sur- 
face pressure term, the expression for virial equilibrium 
can be written as 



3M 



^^1/2 b n , 

81^ + 3^^ 



3a GAP 
1 W 



(7) 



where a and b are two factors of order unity that take 
into account the cloud geometry and the topology of the 
magnetic field. The equilibrium value of is hence given 

by 



B' = 



9aGM 



hb m 



M - 



8 In 2 aC 



(8) 



Assuming a=1.17 and 6=0.3 (McKee et al. 1993) and pos- 
ing M — A/cD we can compute the values of B required 
to support our clumps: these are a few mC Such values 
are comparable to those obtained by Lai et al. (2001) who 
measured the magnetic field in molecular clouds with den- 
sities similar to those of our clumps using interferometric 
polarisation maps. 

We thus propose a scenario in which the clumps are 
marginally stable due to the contribution of the mag- 
netic field. In this case the clump will contract on the 
timescale ^ad determined by ambipolar diffusion. This 
can be related to the free-fall time, iff, by means of 



Eq. (51) of McKee et al. (1993): ^ad 8.5 iff, thus ob- 
taining tAD(yr) = 2.9 10* [nH2(cin~^)]~°'^- The mass ac- 
cretion rate onto the embedded stars can be hence com- 
puted from the ratio Mace = -AfcD/iAD and plotted as 
a function of the source luminosity, as done in Fig. ^ in 
this figure we have also reported the estimates obtained 
from the C^^O data of Hofner et al. (2000). The accretion 
rate increases with source luminosity: this is consistent 
with the finding that the star formation rate in a cloud 
is roughly proportional to the bolometric luminosity ac- 
cording to the expression M*(M0yr~-'^) ~ 4 IO^^'^L(Lq) 
(Plume et al. 1997). Therefore, it is reasonable to expect 
that also Mace satisfies a similar relationship. We hence 
fitted the data in Fig. ^ assuming Mace oc Lfir, thus 
obtaining Mace/ipiR = 2.1 10"* M^yv-'^L'^. The ra- 
tio Af,/Maee — 1-9% is the star formation efficiency in 
the clumps. On the other hand, a lower limit to the same 
value can be obtained from the ratio between the mass 
of the early-type star ionising the embedded UC Hii re- 
gion and the mass of the corresponding clump: this is ap- 
proximately 30 M0/5OOO Mq ~ 0.6%, consistent with the 
previous estimate. It is also worth noting that the mass 
accretion rates are of the order of '^10"^ A/0yr"^: val- 
ues such high could support theories which predict that 
high-mass stars can form through accretion (Behrend & 
Maeder 2001). 

In conclusion, although other explanations are possi- 
ble, we cannot exclude the possibility that the molecu- 
lar clumps traced by the CII3C2II transitions are on the 
edge of gravitational collapse. In support of this hypothe- 
sis we note that evidence for infall in some of our clumps 
has been found by other authors on the basis of high an- 
gular resolution observations in various molecular tracers 
(Hofner et al. 1999; Maxia et al. 2001). In Sect. [4.2. 2| we 
will demonstrate that also the density structure of the 
clumps favours the scenario presented above. 



4.2. Internal structure of the clumps 

We have found that the clumps seen in the CH3C2H lines 
have mean diameters of 0.9 pc, kinetic temperatures of 
~40 K and H2 volume densities of 10^-10^ cm~'^. It is in- 
teresting to compare these values with those derived from 
other tracers. 

For the same quantities, Hofner et al. (2000) find 
slightly different mean values from their C^^O observa- 
tions, namely diameters of '^1.1 pc, kinetic temperatures 
of ^^20 K and H2 volume densities of '^10'* cm"'^. This 
means that the C^^O emission arises from larger, colder, 
and less dense regions than CH3C2H. This is not surpris- 
ing as the C^^O is an order of magnitude more abundant 
than CH3C2H and hence likely to trace lower density gas. 

On the other hand, the C'^^S observations of Cesaroni 
et al. (1991) and Olmi & Cesaroni (1999) provide greater 
densities and smaller diameters than our CH3C2H mCcL- 
surements, namely ^10^ cm~'^ and ~0.4 pc. 
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Source 


D 




-'"CD 


"■H2 




(pc) 


(A/q) 


(Afo) 


(lO^cm^^) 


G5.89 


0.6±0.2 


1800±300 


4500±1200 


8±2 


G9.62 


0.8±0.2 


1000±400 


5000±1000 


5±1 


G10.47 


0.6±0.1 


2700±1000 


6000±2000 


10±3 


G10.30 


0.8 


1300±400 


2000±1000 


20±10 


G10.62 


0.8±0.2 


3000±1000 


8000±2000 


2.5±0.8 


G19.61 


0.5 


1800±500 


1500±500 


7±2 


G29.96 


1.6±0.3 


1300±400 


7000±2000 


0.7±0.2 


G31.41 


1.2±0.3 


2500±600 


5000±2000 


3±1 


G34.26 


0.8±0.2 


2400±700 


6000±1500 


6±2 


G45.47 


1.4±0.3 


2200±800 


6000±2000 


0.9±0.3 


W51D 


0.9 


1700±600 


7000±3000 


50±20 


IRAS 20126 


0.18±0.06 


90±30 


75±20 


5±1 



calculated assuming a CH3C2H abundance of 2 10 



Finally, Kurtz et al. (2000) compiled a list of the prop- 
erties of hot molecular cores (see their Table 1): such cores 
are present in most of our sources and represent the inner- 
most, densest parts of the clumps seen in CH3C2H, with 
typical temperatures of ~100 K and diameters ^0.1 pc. 

In the following we make use of all these information 
to derive the distribution of the physical parameters in the 
clumps. 

4.2.1. Temperature structure 

If we assume that the gas of the clump is heated by an 
embedded high-mass star, the gas temperature, T, at a 
distance R — 9 d from the star is given by the expression 
(see e.g. Plume et al. 1997 and references therein) 



(9) 



with d distance to the source, L luminosity of the embed- 
ded star, and 9 angular distance from the star. In order 
to verify whether such a relation holds for our clumps 

and to compute the values of C and q, in Fig. we plot 
1/2 

Log(T) versus Log[Lp\^/ {9 d)]: note that both quantities 
are distance independent. The distribution can be fitted 
with a straight line corresponding to q — —0.54 ±0.06 and 
C — 10 ± 2, when T is expressed in K, Lfir in Lq, d in 
kpc, and 9 in arcsec. 

Equation (j^ implies that for a given source T cx 0'. 
Therefore, one can also derive q from the ratio between the 
temperatures measured in two different tracers and the 
ratio between the corresponding diameters, which satisfy 
the relation 




Ty 



i?x 



(10) 



with X,Y indicating that the value has been obtained from 
CH3C2H, C^^O (Hofner et al. 2000), or taken from Tabic 1 
of Kurtz et al. (2000). Under the assumption that q is 
exactly the same for all sources, Eq. (|l^ is totally source 



ao 
o 



0.5 1 1.5 2 

T'Og[T''4c (T;/^ i)/0(arcscc)] 

Fig. 8. Plot of the kinetic temperature versus the quantity 
Likpc/0, where Likpc is the source luminosity scaled to the 
distance of 1 kpc and 9 the angular radius at which the 
temperature has been measured. Filled circles and squares 
indicate respectively estimates based on CII3C2H (this 
study) and C^^O (Hofner et al. 2000), while empty tri- 
angles correspond to the values of Table 1 of Kurtz et al. 
(2000). The straight line represents a least square fit to 
the data, with slope 0.54 ± 0.06 and intercepta 1 ± 0.09 



independent and one can plot all possible ratios for all 
sources in the same figure. The advantage of this approach 
with respect to that previously used is that we do not need 
any luminosity estimate. The result is shown in Fig. ^: the 
slope in this case is q = —0.56 ± 0.06, consistent with the 
previous estimate. 
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o 



-0.5 



Log(0y0,) 

Fig. 9. Plot of the temperature ratio Tx/Ty versus the 
angular diameter ratio 6x/0y, where X and Y indi- 
cate two estimates among those derived from CH3C2H, 
C"0 (Hofner et al. 2000), and those of Table 1 of Kurtz 
et al. (2000). Circles correspond to ratios between the 
temperatures quoted by Kurtz et al. and those obtained 
from CH3C2H, squares to ratios between estimates from 
Kurtz et al. and from C^^O, and triangles to ratios be- 
tween estimates from CH3C2H and C^^O. The straight 
line represents a least square fit of the type T oc 0*, with 
q — —0.56 ± 0.06. The correlation coefficient is -0.93 

We conclude that the temperature in the clumps is ap- 
proximately described by the law T oc i?^"-^^. Such a law 
is ste eper t han expected on the basis of the discussion of 
Sect. 3.4.2, where we conclude that no temperature gra- 
dient is seen in the gas traced by the CH3C2H emission. 
However, in order to reveal this gradient we had to com- 
pare the temperature measured on the large scale (in the 
CH3C2H and C^^O emission) with that on a very small 
scale (in CH3CN). This means that while T oc R-°-^^ is a 
reasonable description of the "global" temperature varia- 
tion from ^0.1 pc to ~1 pc, it cannot be applied rigorously 
at all radii in the clump. In particular, the temperature 
profile is likely to be shallow in the outer region (for the 
reason discussed in Sect. 3.4.2) and steeper close to the 
center. 

Note also that q=-0.55 is significantly greater than the 
slope expected for an optically thin clump heated by an 
embedded energy source. In fact, in this case for a dust 
absorption coefficient (xiy^ one finds T oc the 
steepest slope is obtained for /3 = 1 and is equal to -0.4. 
A plausible explanation for q being less than -0.4 is that 
part of the clump is optically thick. Indeed, in some of 
our objects, on scales 10 times smaller than the clumps, 
Cesaroni et al. (1998) found optically thick, hot cores with 



temperature profiles T oc R^^ ''^ . These cores are just the 
regions which the temperature estimates by Kurtz et al. 
(2000) used in our plots refer to. We thus conclude that 
one may read the slope derived by us as a "compromise" 
between a flatter temperature profile in the optically thin 
clump - traced by C^''0 and CII3C2II - and the steeper 
profile (T oc R'^-'^^) in the hot core - measured e.g. in 
NH3 or CH3CN. 



4.2.2. Density structure 

For star forming clumps, theoretical models predict den- 
sity profiles of the type n oc R^ . In our case, it is diffi- 
cult to test this prediction source by source, as only few 
measurements for each source at different radii are avail- 
able, namely those from CH3C2H, C^^S, C"0, and the 
density estimates obtained from Table 1 of Kurtz et al. 
(2000). However, the ratio between the H2 densities mea- 
sured with different tracers is equal to 



(11) 



with nx,Y density measured with tracer X,Y=CH3C2H, 
C^^O, C'^^S, and mm contimmm (from Kurtz et al. 2000), 
and i?x.Y radius of the spherical region traced by X,Y. 
One can demonstrate that if Eq. ( pT| ) holds for the volume 
density at radius i?, then it must hold also for the mean 
density inside R, and vice versa. This result is important 
because it makes possible to use the mean densities given 
in Table ^ to derive an estimate of p. 

As already done for the temperature profile, we assume 
that p is the same for all sources. Under this hypothesis, 
one can plot Log(rix/'^Y) versus Log(9x/0Y) for all pos- 
sible ratios and all sources: this is shown in Fig. The 
correlation is excellent (the correlation coefficient is -0.97) 
and well reproduced by a linear fit with p = —2.59 ± 0.07. 
This result is consistent with the density profiles predicted 
by star formation models (Shu et al. 19 87; Li 1999) and 
supports the idea discussed in Sect. [4.l| that our clumps 
are marginally stable. Noticeably, in W3(H20) - an ob- 
ject of the same type as those observed by us - an even 
steeper density law has been derived by Mauersberger et 
al. (1988): they find indirect evidence for p ranging be- 
tween -3 and -2.8. This is in reasonable agreement with 
our finding. 

It is worth noting that the previous density law applies 
to both the core (empty points in Fig. Q and the sur- 
rounding clump (filled points) independently. This result 
is not entirely consistent with the findings of Hatchell et 
al. (2000), who fit their continuum maps of similar clumps 
around UC Hii regions with shallower density profiles. 
However, their fits are model dependent, whereas our re- 
sults have been obtained under relatively simple assump- 
tions, which in our opinion make them more reliable. 

On the other hand, we know that massive stars have al- 
ready formed inside the clumps, which means that gravita- 
tional collapse is going on in the central regions. According 
to Shu et al. (1987) this occurs within a radius i?2(pc) = 
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Fig. 10. Same as Fig. || for the H2 density. Filled sym- 
bols: circles correspond to ratios between densities from 
C^^O and C^^S, squares to ratios between densities from 
C^'^S and CH3C2H, and triangles to ratios between den- 
sities from C^^O and CH3C2H. Empty symbols: circles, 
triangles, and squares correspond respectively to ratios 
between densities estimated from Table 1 of Kurtz et al. 
(2000) and those derived from C"0, CH3C2H, and C^^S. 
The straight line represents a least square fit of the type 
n oc 9^, with p — —2.59 ± 0.07. The correlation coefficient 
is -0.97 



1.72 10-3 M4M0) [Awi/2(kms-i)]-2, where M* is the 
mass of the star created by the collapse: inside this ra- 
dius the density should be (xR~^/^. In our case, the max- 
imum value of i?2 is obtained for ~ 60 Mq and 
At)i/2 — 2 km s"^, which give R2 — 0.026 pc, much less 
than the minimum radius measured in any of the tracers 
used in Fig. |l^: this confirms that the density profile mea- 
sured with such tracers can be ocR~^ , close to our result. 

It is worth noting that n cx R^'^ holds for an isothermal 
sphere, which seems to contradict the result T oc R~^'^^ 



found in Sect. 4.2.1. However, the density gradient is es- 
tablished on dynamical timescales, much longer than the 
radiative timescale needed for the newly born high-mass 
stars to heat up the clumps. As a result, the density pro- 
file is a "remnant" of the primeval structure of the clump, 
whereas the temperature profile refiects the present heat- 
ing by the embedded stars. 

4.2.3. Velocity structure 

Given the low angular resolution of our observations, the 
only possibility to investigate the velocity field in the 
clumps is to study the variation of the line shape as a 



function of the corresponding excitation energy: transi- 
tions from higher energy levels are expected to arise from 
inner regions since the temperature increases towards the 
center of the clumps, as previously discussed. This makes 
possible to sample the clumps at different radii by means 
of the different K lines of CH3C2H. 

To this purpose, we have relaxed the assumption of 
equal line FWHM and peak velocity used in the line fits 
illustrated in Sect. |2.3| and we have fitted each K compo- 
nent independently from the others. Then we have consid- 
ered the variation of the line FWHM and peak velocity as 
a function of excitation energy. To this purpose, for each 
source we have computed the ratio between the FWHM 
of each line and that of a reference line, and the difference 
between the velocity of each line and that of the same 
reference line. The latter was the (6-5) K—2 transition, 
which is strong and does not overlap with other K lines. 
Such ratios and differences have been plotted as a func- 
tion of the ratio between the energies of the correspond- 
ing transitions and that of the reference line: the result is 
shown in Fig. 11. Clearly, no trend is seen for the velocity: 
this suggests that all CH3C2H lines trace the whole emit- 
ting region, as expected for optically thin transitions. On 
the contrary, the FWHM ratio shows a slight tendency to 
increase with the energy of the transition. 

This fact cannot be explained with thermal broaden- 
ing of the lines, because temperature changes at most from 
~10 K on the clump surface to ^100 K in the inner re- 
gions: the corresponding contribution to the line width 
is respectively 0.1 km s^^ and 0.34 km s^^, which for a 
non-thermal line width of at least '^3 km s^^ determines a 
variation of 0.6% of the line width. The latter is negligible 
to our purposes. 

Another explanation for the observed line broadening 
at high energy might be an increase of optical depth for the 
high-energy transitions. It is possible to demonstrate that 
for an optical depth t at the line centre, the ratio between 
the line FWHM, Aw, and the corresponding value in the 
optically thin limit, Awq (i.e. the intrinsic line width), is 
given by the expression 



At; 
A^ 



In 2 



In 



ln2-ln(l + c-^) 



(12) 



From this, one can see that r ~ 4 is needed to ob- 
tain Au/Aup — 1.5, i.e. the maximum ratio seen in the 
top panel of F ig, pi] . As already discussed at length in 
Sects. |3.2| and 3.4.1, such a large value must refer to a 
very small region, which cannot contribute to the CH3C2H 
emission observed by us on a much larger scale. Therefore, 
line broadening is unlikely to be due to optical depth ef- 
fects. 

In conclusion, we believe that the most plausible expla- 
nation for Fig. |l^ is that the velocity dispersion increases 
with energy of the transition and hence towards the centre 
of the clumps. Various mechanisms can explain this trend: 
keplerian rotation, infall motions, and turbulence caused 
by winds/outflows from the embedded, young high-mass 
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Fig. 11. Top panel: plot of the ratio between the FWHM 
of a CH3C2H K line and the FWHM of the (6 - 5) K^2 
line, versus the corresponding excitation energy ratio. The 
line connects the mean values of the ratios obtained for 
each transition: one can see a slight increase of line width 
with energy. Bottom panel: same as top panel, for the 
difference in peak velocity between each K line and the 
(6 - 5) A' = 2 line. Unlike the FHWM, no trend is present 
and the peak velocity is approximately the same for all 
transitions 



— we have estimated the diameters of the CH3C2H emit- 
ting regions by comparing the main beam brightness 
temperatures measured with the 30-m and the FCRAO 
telescopes: a typical diameter turns out to be ~ 0.9 
pc, of the same order of those obtained with other 
medium density tracers (Cesaroni et al. 1991; Hofner 
et al. 2000); 

— we have assessed that the CH3C2H molecule is in LTE 
and derived an accurate estimate of the temperature of 
the molecular clumps associated with these high-mass 
star forming regions: the temperatures range from 30 
to 56 K; 

— we have derived CH3C2H column densities of 10^^- 
10^^ cm~^ corresponding to H2 volume densities of 
~10^-106 cm~3; 

— by comparing our results with those obtained by other 
authors using different tracers (Cesaroni et al. 1991; 
Hofner et al. 2000; Kurtz et al. 2000) we find that tem- 
perature and density depend on the distance from the 
clump centre, R, according to the relations n oc R~^-^ 
and T oc i?^*^'^: these are consistent with marginally 
stable clumps heated by the embedded massive stars; 

— the masses of the clumps turn out to be ~3 times 
greater than the corresponding virial estimates ob- 
tained ignoring magnetic fields and external pressure: 
we suggest that this might indicate that the clumps are 
on the verge of gravitational collapse and that mag- 
netic fields might play an important role in stabilising 
the clumps; we find that the magnetic fields needed for 
virial equilibrium are of the order of 1 mG; 

— the dynamical timescales on which the clumps may 
collapse are of ^ 10^ yr: these imply mass accretion 
rates of - 10"^ Mgyr'^ 

— the line width shows a tendency to increase for in- 
creasing excitation energy, and hence for decreasing 
distance from the centre: a plausible explanation is 
that such a line broadening is due to collapse of the 
inner regions. 
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stars. Although the present observations cannot discrimi- 
nate among these, it is worth noting that the infall hypoth- 
esis is con sistent with the scenario proposed in Sects. 4T 
and |4.2.2|. 



5. Conclusions 

We have used the IRAM 30-m and FCRAO 14-m tele- 
scopes to observe a sample of 12 UC Hii regions in the 
J=6-5, 8-7, and 13-12 transitions of CH3C2H. Emission 
has been detected in all of the objects and the following 
results have been obtained: 



References 

Behrend, R. & Maeder, A. 2001, A&A 373, 190 

Bergin, A.E., Goldsmith, P.P., SneU, R.L. & Ungerechts, H. 

1994, ApJ, 431, 674 
Cesaroni, R., Hofner, P., Walmsley, C.M. & Churchwell, E. 

1998, A&A, 331, 709 
Cesaroni, R., Walmsley, CM., Kompe, C. & Churchwell, E. 

1991, A&A, 252, 278 
Cesaroni, R., PeUi, M., Jenness, T. et al. 1999, A&A, 345, 949 
Chini, R., Kreysa, E., Mezger, P.G. & Gemiind, H.-P. 1986, 

A&A, 154, L8 

Churchwell, E., Walmsley, C.M. & Cesaroni, R. 1990, A&AS, 
83, 119 



Fontani et al.: CH3C2H in UC Hll regions 



15 



DubruUe, A., Boucher, D., Burie, J. & Demaison, J. 1978, J 

Mol Spec, 72, 158 
Erickson, N.R., Goldsmith, P.F., Novak, G. ct al. 1992, IEEE 

Transactions on Microwave Theory and Techniques, 40, 

N°l 

Hatchell, J., Puller, G.A., Millar, T.J. et al. 2000, A&A, 357, 
637 

Hofncr, P., Kurtz, S., Churchwcll, E., Walmsley, CM. & 

Cesaroni, R. 1996, ApJ, 460, 359 
Hofner, P., Wyrowski, P., Walmsley, CM. & Churchwell, E. 

2000, ApJ, 536, 393 
Hofncr, P., Peterson, S. & Cesaroni, R. 2000, ApJ 514, 899 
HoUis, J.M. 1982, ApJ, 260, 159 

Hunter, T.R., Churchwell, E., Watson, C. et al. 2000, AJ, 119, 
271 

Irvine, W.M., Goldsmith, P.P. & Hjalmarson, H. 1987, 
Interstellar Processes, cd. HoUcnbach D.J., Tlironson H.A., 
561 

Kurtz, S., Cesaroni, R., Churchwell, E., Hofner, P. & Walmsley, 
CM. 2000, Protostars and Planets IV, ed. Mannings V., 
Boss A., Russol S. (Tucson: Univ. of Arizona Press) 
Kutner, M.L. & Ulich, B.L. 1981, ApJ, 250, 341 
Lai S., Crutcher R.M., Girart J.M.& Rao R. 2001, ApJ, 561, 
864 

Li, Z. 1999, ApJ, 526, 806 

MacLaren, I., Richardson, K.M. & Wolfendale, A.W. 1988, 
ApJ, 333, 821 

Mauersberger, R., Wilson, T.L. & Henkel, C. 1988, A&A 201, 

123 

Maxia, C, Tcsti, L., Cesaroni, R. & Walmsley, CM. 2001, 
A&A 371, 287 

McKee, C.F., Zweibel, E.G., Goodman, A.A. & Heiles, C. 1993, 
in Protostars and Planets HI, eds. E.H. Levy & J.I. Lunine, 
p. 327 

Mczgcr, P.G., Zylka, R. & Wink, J.E. 1990, A&A, 228, 35 
Olmi, L., Cesaroni, R. & Walmsley, CM. 1993, A&A, 276, 489 
Olmi, L. & Cesaroni, R. 1999, A&A, 352, 266 
Plume, R., Jaffe, D.T., Evans II, N.J., Martin-Pintado, J. & 

Gomez-Gonzales, J., 1997, ApJ, 476, 730 
Shu, F.H., Adams, F.C & Lizano, S. 1987, ARA&A, 25, 23 
Spitzer, L., 1978, Physical processes in the interstellar medium, 

ed. John Wiley & sons. Inc. 
Townes, CH. & Schawlow, A.L. 1975, Microwave spectroscopy, 

ed. Dover, Now York 
Wilson, T.L. & Rood, R. 1994, ARA&A, 32, 191 
Wood, D.O.S. & Churchwcll, E. 1989a, ApJS, 69, 831 
Wood, D.O.S. & Churchwell, E. 1989b, ApJ, 340, 265 



